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, Abstract 

We describe solutions of the Klein-Gordon equation which are spherically symmet- 
ric and localized, and may be regarded as massive particles without charge or spin. 
' ^ | The proposed model, which is based on the action for a complex scalar field minimally 

^ I ■ coupled to the electromagnetic and gravitational fields, contains no adjustable param- 

eters and predicts five particle species with masses of the order of the Planck mass. 
f"^ , These particles appear to be candidates for dark matter. 

\o ■ 

in ! 1 Introduction 

O ■ 

| A number of models describing particles as solitons — stable self-bound concentrations of 

(3^. field energy — have been proposed. In an early paper Rosen |14| obtained soliton solutions 

| from the interaction of a complex scalar field and the electromagnetic field. Later Cooper- 

bJQ' stock and Rosen [S] obtained soliton solutions by coupling a complex scalar field to both 

electromagnetism and gravity. In the meantime Wheeler |16| . using an entirely different 
•'■h . theoretical framework, obtained soliton solutions (which he dubbed geons) from a purely 

/\ ' classical model of electromagnetism coupled to gravity. More recently Moroz, Penrose and 

Tod ^21 obtained soliton solutions of the Schrddinger-Newton equations. See El 12] 

and references therein for a window into the literature. 

The model presented here is similar to that of Cooperstock and Rosen. 1 What differs is 
the treatment of boundary conditions and the introduction of certain constraints. Cooper- 
stock and Rosen require all fields to be finite and continuous at the center of the particle and 
they integrate the field equations outward. Parameters are adjusted until the asymptotic 
wave function vanishes. The resulting model contains several adjustable parameters. 

In contrast, we impose asymptotic boundary conditions and integrate the field equations 
inward, toward the center of the particle. And we require the asymptotically measured mass 
and charge to equal the mass and charge parameters of the action. The resulting model 
contains no adjustable parameters. 

In Section [21 we write down the action for a complex scalar field minimally coupled to 
the electromagnetic and gravitational fields, and use the principle of stationary action to 
derive the corresponding field equations (the Klein-Gordon, Einstein and Maxwell equa- 
tions). Then we propose a stationary spherically symmetric trial solution for the metric 
tensor, wave function and vector potential and substitute it into the field equations. We 
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1 Oxii Eqs. fl^-flS} are equivalent to Eqs. (5.3)-(5.6) of 0. 
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obtain a system of four coupled nonlinear differential equations, as well as two integral con- 
straints related to the conservation of charge and energy. In Section|21we search for localized 
particle-like solutions. After calculating the asymptotic behavior of the equations and im- 
posing asymptotic boundary conditions, we are left with a single adjustable parameter — the 
system electric charge. By again appealing to the principle of stationary action we find that 
the charge must vanish. Hence, we arrive at an eigensystem with no adjustable parameters. 
Using numerical methods we calculate the eigenmodes and find five massive particle species 
without charge or spin. We name these particles quantum geons and show that the space- 
time curvature diverges at the center of each geon. In Section0|we explore the implications 
of this singularity and argue that it is benign. In Section [5] we discuss our results. 



2 Field equations 

Consider the action 

S = ^r I [ R - flM,f ^ + iC*^*:/* + * : ^) - ^ 2 **] Vd 4 x (i) 

I" 7 '" Jv 

for a complex scalar field \& within a four-dimensional volume V with spacetime curvature 
R, where a colon denotes the generalized covariant derivative p- 167] 

¥.„ = + iQA^ , (2) 

a semicolon denotes the covariant derivative of general relativity, an overline denotes complex 
conjugation, Q is electric charge, is the electromagnetic vector potential, 

F^tx = Aj_i : u (3) 

is the electromagnetic field tensor and M is mass (assumed positive). All quantities are 
expressed in natural units where the speed of light, Planck's reduced constant (K), Newton's 
gravitational constant and Coulomb's electrostatic constant are unity. Our notation and 
conventions are detailed in Appendix lAl 

Make small arbitrary variations SS in the action by making small arbitrary variations 
Sga/3, <5A^, S^ and in the fields, and neglect boundary terms arising from integration by 
parts. Then impose the condition SS — for stationary action. The terms proportional to 
5^ and 6*f? yield the Klein-Gordon equation 

+ M 2 * = (4) 

and its complex conjugate. The terms proportional to Sg a p yield the Einstein equations 

G af) = -87rT Q/3 , (5) 

where 

G a/3 _ R a/3_ l g a/3 R (g) 

is the Einstein tensor and 

16-7T V / 
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is the energy-momentum tensor. The terms proportional to <5A M yield the inhomogeneous 
Maxwell equations 

F"".„ = 4tt_P , (8) 

where 



16tt v ; 
is the electromagnetic current vector. The homogeneous Maxwell equations 



(9) 



(10) 



follow immediately from the definition of F „„. 

Let the volume V in Eq. Q include all of space over some arbitrary time interval T. 
We will later impose boundary conditions which require R, A p and \I/ to vanish asymptoti- 
cally. Thus we are justified in neglecting boundary terms when integrating by parts in the 
expression for SS. 

When Eqs. (@J, (0 and are satisfied the action is stationary and the scalar curvature 

is 

R = 8ttT = 2Af 2 ** - vP 7 *^ , (11) 
where T = T M M . Substituting for R in Eq. gives 



as the stationary value of the action. 



(12) 



2.1 Trial solution 

Our goal is to solve Eqs. (0, 10 and for the fields g M „, '5 and A M . In order to make the 
calculations tractable we consider a static spherically symmetric spacetime with spherical 



polar coordinates x 



t, x l 



r, x 



9 and x 



and metric tensor 



( u 








° ^ 





— V 
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V o 








-r 2 sin 2 6> / 



(13) 



where u and v are real functions of r. The corresponding spherically symmetric wave 
function and vector potential are assumed to have the forms 

-iut 



f = Re 



(14) 



and 



QA M = 



(15) 






V o ) 

where R and $ are real functions of r and to is a real constant. 

Expressions (|13|) . i|14|) and (|15fl for g p „, W and A M comprise a trial solution of field 
equations 0, and 10 and allow us to calculate all quantities of interest in terms of R, 
u, v and <&. We will show that the field equations reduce to a set of four coupled nonlinear 
differential equations in R, u, v and Then we will solve these equations subject to 
physically plausible boundary conditions and constraints. The fact that the field equations 
are, indeed, satisfied by the assumed trial solution is its ultimate justification. 
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2.2 Einstein, energy-momentum and electromagnetic tensors 

Given trial solution I|14l) and (|15fl we can calculate the quantities in field equations 

(@J, JSJ and ©. The nonzero elements of the Einstein and energy-momentum tensors are 



Gi 1 



G 2 2 



v' 


1 


1 






r 2 v 


u' 


1 


1 


ruv 




r 2 v 


u" 


{u'f 


u'v 


2uv 


Au 2 v 




G 2 2 







and 



V = 
T 2 2 = 



T 3 - 
I 3 — 



Af 2 i? 2 


+ (« 


-$) 2 i? 2 


(i?') 2 


16tt 




167TM 


167TW 


M 2 i? 2 


( W 


-$) 2 i? 2 


(i?') 2 


16tt 




167TM 


167TU 


M 2 i? 2 


( W 


-$) 2 i? 2 


(flf 


16tt 




167TM 


167TW 



2ruu 2rv 2 



(*') 



(16) 



8irQ 2 uv 

m 2 

8irQ 2 uv 

m 2 

8irQ 2 uv 



(17) 



where primes denote differentiation with respect to r. The contraction of the energy- 
momentum tensor is 



T = 



M 2 R 2 (w-$) 2 i? 2 {R'f 



4-7T 8iru 8irv 

The nonzero elements of the electromagnetic field tensor are 

QFqi = -QFio = $' 

and the only nonzero element of the electromagnetic current vector is 



J° = 



Q (u - $) R 2 



8nu 



We also note that 



Q 2 F"%„ = - 



2(*') 



(18) 
(19) 

(20) 
(21) 



2.3 Differential equations 

The Klein-Gordon equation (0} becomes 



u 
2ii 



v 

2^ 



(uj - $) 



M 



= 0. 



(22) 



There are four Einstein equations JSJ), one for each nonzero component of . However, 
the equations corresponding to G2 2 and G3 3 are identical. Furthermore, they can be derived 
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by combining the Go°, Gi 1 and Klein-Gordon equations. So we retain the Go°, Gi 1 and 
Klein-Gordon equations, and discard the G2 2 and G3 3 equations. It will prove convenient 
to work with the equations obtained by adding and subtracting the Go° and Gi 1 equations: 



and 



u' v' 1 — v 
2u 2v r 



u v rv 
2u + 2v ~ Y 



M 2 R 2 



2 Q 2 uv 
(w-$) 2 i? 2 (i?') 2 



U V 

The inhomogeneous Maxwell equations JHJ yield the Poisson equation 

Q 2 R 2 v 



2 v! v'\ 
] $' 

r 2u 2v , 



(U3 - = . 



(23) 



(24) 



(25) 



The four coupled nonlinear differential equations \\T1\\ through l|25|l are to be solved for 
the functions R, u, v and $ subject to physical constraints and boundary conditions. 

2.4 Conserved quantities 

In this section we derive expressions for the electric charge, energy, angular momentum and 
stationary action of our model. 



2.4.1 Charge 

The electromagnetic current vector satisfies the conservation law 

-i% = (JV),, = o. 

The corresponding conserved quantity 



l°^drdQd<t> 



(26) 



(27) 



is the electric charge, where k q is a real constant (to be determined later). Substitute J° 
from Eq. J5UJ) and let 

V = ^ 2 v / H si n6 l (28) 

then integrate over 9 and tj> to get 



I = hi 

2 



(w - *) R 2 



I dr. 



(29) 



Another expression for Q can be obtained as follows. Substituting from Eq. |(HJ into Eq. 
(|T7Jl and using the identity 

P"VV=(RVL (30) 



gives 



(31) 
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For our static spherically symmetric metric the integrand is 



(F°V) 



dr V Quv 



V 



•,2 _ kq 



47r ./„ dr 



V )dr d9 < 



(32) 
(33) 



Substituting expression (|28[) for then integrating over 9 and </> gives 

Q 2 = 



o 

2 (oc) : 



dr 



where 



! (r) 



(34) 



(35) 



and we assume (to be justified later) Q 2 (0) = 0. The quantity Q 2 (rb) — Q 2 (r a ) is the square 
of the charge contained in the region between r a and n,. 

2.4.2 Energy 

In our static spherically symmetric spacetime the vector field 



satisfies the Killing equation 
If we let 



( 1 \ 




V o / 



J £ — S \u 



then the anti-symmetry of implies 
and 



= 0. 



(36) 

(37) 
(38) 

(39) 
(40) 



Thus is a conserved current associated with the time invariance of the metric. The 
corresponding conserved quantity 



J CO 



(41) 



G 



is the system mass (energy), where is a real constant (to be determined later). Calculation 
of from Eqs. and l|38|) for our static spherically symmetric spacetime gives 



/Ro°\ 




V o J 



But 

so Eq. (|41|) becomes 



-8tt (To 



o 1 



T 



M = -87rfc ? I (T ° - ±T) V dr d6 d(j> . 

J OO 



(42) 



(43) 



(44) 



Substituting expressions fTT)) . (fH^I and (|2"%f for To°, T and yj into (|T4^I . then integrating 
over 9 and <p gives 



M = -47rfc ? 



(w - $) 2 i? 2 M 2 i? 2 ($') 2 



(45) 



Another expression for M can be obtained as follows. Substituting from Eq. (|39(l into 
Eq. iglj gives 

(46) 



M = ki f (£ 0; V),„ drd0( 

t/ OO 



For our static spherically symmetric metric the integrand is 



d 



Substituting expression (|28|l for then integrating over 9 and <fi gives 



(47) 
(48) 



M = 



o dr 
M(oo) , 



where 



M(r) 



_ 27r/c 5 r 2 u'^/|uw| 



(49) 



(50) 



and we assume (to be justified later) M(0) = 0. The quantity M(rt) — M(r a ) is the mass 
contained in the region between r a and r\,. 
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2.4.3 Angular momentum 

Based on the Killing field 






V i J 



(51) 



and using an approach similar to that of Section \'2. 4 . 21 one can show that the system angular 
momentum is zero, as expected for a spherically symmetric model. 



2.4.4 Action 

From Eqs. lfT2)l. (fl4^ . lt2"T|l and l(2"%|l the stationary value of the action is 



T 

S =2 



APR 2 



(gr 

Q 2 uv 



(52) 



3 Solution of the field equations 
3.1 Solution at infinity 

Our goal in this section is to determine the asymptotic (r — > oo) expressions for the fields R, 
u, v and <&. We seek a particle-like solution with charge and energy distributions localized 
near the origin and we impose the following boundary conditions: 

1. The wave function \E' and all components of the vector potential A M asymptotically 
vanish. 

2. Spacetime is asymptotically flat and local measurements by an asymptotic observer 
agree with the predictions of flat spacetime physics. Specifically, the asymptotic be- 
haviors of u and v reproduce Newtonian gravity, and the asymptotic expressions for 
i> and R satisfy the flat spacetime Poisson and Klein-Gordon equations. 

These boundary conditions will guide us to candidate expressions for the asymptotic fields. 
We will then verify by direct substitution that these candidates asymptotically satisfy field 
equations i|2"2"jl through 

The asymptotic expressions 



u = l-IMr- 1 + Q 2 r~ 
l 



I! 



U 



(53) 



correspond to the (Reissner-Nordstrom) spacetime structure outside a charged nonrotat- 
ing black hole. They yield flat spacetime (Newtonian) gravity as r — > oo, as required by 
boundary condition |5J They also imply [see Eq. (|50|l ] 



M(oo) = -Ank^M . 



Comparing this with Eq. I|49|l gives 



h = —. 



1 

4tt 



(54) 



(55) 



8 



The asymptotic expression 

$ = Q 2 r- X (56) 

follows from boundary conditions 2] and [21 and the flat spacetime laws of electrostatics. 
Expressions (|B3|) and ij5f)|) imply [see Eq. 1(55)1] 

Q 2 (oo) = fc g Q 2 . (57) 

Comparing this with Eq. 1)34)1 gives 

kq = l. (58) 

The asymptotic form of R follows from boundary condition [21 and the asymptotic flat 
spacetime Klein-Gordon equation [obtained by letting u, v — > 1 in Eq. i)22)l ] 



i?" + - ii' 

r 



(w-$) 2 -Af 2 i? = 0, (59) 



where $ is given by 1)56)) . The solution which asymptotically vanishes (boundary condition 
EJ is 

^=%^> (60) 



where i?^ is a real constant, 



and 



° = —r- » ( 61 ) 



fe = VM 2 - lu 2 (62) 

with w 2 < M 2 . 

We now have asymptotic expressions for the helds which satisfy the boundary conditions. 
But we still need to verify that these expressions are consistent with our field equations. 
This is accomplished by substituting the asymptotic expressions for R, u, v and <& into Eqs. 
(|22|l through (|25|l and verifying that the dominant terms satisfy the equations as r — ► oo. 
This procedure reveals that Eqs. 1)23)) through (|25ll are, indeed, asymptotically satisfied. 
However, Eq. 1)221) is satisfied only if we let (see Appendix lB|l 

Q2U M ' M 2 - 2. 2 ) , (63) 



k k 

which differs from l|61|l . In order to satisfy boundary condition [21 we must make 1)61)) and 
()63)1 consistent, which requires 

lu = ±M/V2 (64) 

and [from 1)62)1] 

k = M/V2. (65) 

In summary, based on the imposed boundary conditions we expect the asymptotic solu- 
tion to have the form 

R = R^r-^&e-^/^ 
u = 1 - 2MT- 1 + Q 2 r- 2 
v = u~ 1 

$ = QV" 1 . (66) 
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3.2 Physical constraints 



In addition to satisfying field equations H22|) through 1251) the solution must also satisfy the 
charge and energy constraints of Section 12.41 Using results from Section 13.11 in Eqs. (|29|) 
and 145fl we obtain the constraints 



1 = 




(67) 



and 



1 = — 

M 



± 



M 
V2 



M 2 R 2 



m 2 

Q 2 uv 



(68) 



3.3 Solution strategy 

If we are given values of the three parameters Q, M and Roo and a choice of sign in Eq. (|64(l 
then the asymptotic fields are fully defined by Eqs. (|66J) . A complete solution can be obtained 
by integrating field equations l|22|l through (|25|) from the asymptotic limit (r — > oo) down 
to the origin (r — * 0). Since the two constraints of Section ET2l must be satisfied, only one 
of the three parameters is independent. It will be convenient to take Q as the independent 
parameter and regard M and R^ as functions of Q. 

Once a solution is found for given Q then the action S can be calculated from Eq. I|52|l. 
We may regard the action as a function of Q. The field equations, constraints and boundary 
conditions actually depend on Q 2 , not Q, so we may write the action as S(Q 2 ). Solutions 
of physical interest have stationary action so they must satisfy 

|| = 2QS'(Q 2 H0. (69) 

3.4 Solution for zero charge 

We will consider solutions with Q = 0. These solutions satisfy Eq. (|69|l so they are of 
physical interest. 

With Q = the asymptotic fields l|66f> become 

r -l e -Afr/V2 



R = 




u = 


1 - 






V = 


u- 1 


$ = 


0. 



(70) 



Field equation (|25|l becomes 



*"+(-- f "fV = (71) 
\r 2u 2v J 

which, given the asymptotic expression for $, yields the trivial solution (valid everywhere) 

$ = . (72) 
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Then field equations (|22|) through (|24[l become 



R" + I 1 + _ ^ ) R> + I jL _ i ) M 2 vR = 0, (73) 



and 



r 2u 2i>, 

u' u' 1 - « _ M 2 rvR 2 

2u ~ 2v + r ~ 2 

~M 2 R 2 (i?') 2 



u v rv 
2u + 2v ~ T 



2« 



(74) 
(75) 

(76) 



and constraints l|67|) and l|68[) can be written 

foo 2t>2 

A = 2V2tM W\uv\dr = 

Jo u 

and 

pOO 

±2V2-2-Af / r 2 R 2 ^/\u~^\dr = Q. (77) 
Jo 

Since both M and the integral in l|77[) are positive this constraint can only be satisfied if we 
choose the upper sign. So we need only consider the upper signs in Eqs. (|76(l and Ij77(l . and 
[from Eq. uj = 

Given values of M and i?oo, and starting with the asymptotic fields 1|70J) . Eqs. (|73|1 
through (|75|l can be numerically integrated all the way to the origin. The values of M and 
i?oo can be adjusted until constraints (|76|l and i|77|) are satisfied. 

Figure ^ shows the locus in the M-R^ plane which satisfies constraint 1|77[1 . Figure |2 
plots A versus M along that locus. Since A = when constraint (|76[) is satisfied the five zero 
crossings in Fig. [21 correspond to the solutions we seek. These solutions describe uncharged 
spinless massive particles. Because the masses appear spontaneously as eigenvalues of a self- 
gravitating quantum field ( "mass without mass" ) we borrow terminology from Wheeler 
and dub these particles quantum geons. We will index the five geons by the integer n = 
0, 1, 2, 3, 4 in order of increasing mass. 

Various features of the geons are plotted in Figs. El through and numerical character- 
istics are tabulated in Tabled The probability density is 

P = - _ V = NR 2 y/ , (78) 

where the normalization factor 

M 

N = (79) 

8tt( v / 2-1) 

is determined from Eq. 177|) . The radial probability density 

P r = 4:T:Nr 2 R 2 y^\uV\ (80) 

is plotted in Fig.0 The stationary value of the action calculated from Eqs. (|52|) and (|77|l is 

(y/2 - 1) TM , 
S=^ ^ (81) 

and the action per cycle (T = 2tt/u) is (2 — ^/2)n for each of the five geons. 
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Figure 1: Locus in the M-Roo plane which satisfies constraint (|77|l . 
0.8 
0.6 
0.4 

< 

0.2 


-0.2 

0.4 0.6 0.8 1 1.2 1.4 1.6 
M 

Figure 2: Left-hand side of constraint (|76|l versus mass for the locus of points of Fig.Q] The 
five zero crossings correspond to solutions of field equations l|73|) through (|75|l which satisfy 
boundary conditions l|70(l and constraints H76JI and 1|77|L There are no zero crossings in the 
neighborhood of M = 1.25. 

Table 1: Numerical characteristics of the geons. Values of M and i?oo are believed to be 
accurate to within ±1 of the least-significant digit. 



n 


nodes in R 


M 


Roo 








0.605 


1.51 


1 


1 


0.729 


2.22 


2 


1 


0.866 


4.26 


3 


2 


1.000 


9.3 


4 


2 


1.077 


15.6 
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Figure 3: Radial wave function versus radius. Curves labeled by geon index n. 
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Figure 4: Metric function u versus radius. Curves labeled by geon index n. 
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Figure 6: Radial probability density versus radius. Curves labeled by geon index n. 
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3.5 Solution at the origin 

The results of Section 13.41 suggest that in the neighborhood of the origin 



R 

u 
v 



R - 
v 2 r 2 



Re lnr 



(82) 



where Rq, Rg, uq and v 2 are real constants and uq,v 2 > 0. These expressions satisfy Eqs. 
(I73[l through l|75|l as r — > provided 



Re = ±V2. 



(83) 



Now that we know the behavior of the fields as r — > we are in a position to justify two 
earlier assumptions. From expression l|35|) for Q 2 (r) we obtain 



lim Q\r) 

r^Q 



lim 



0. 



justifying the assumption Q 2 (0) — 0. And from expression H5()|) for M(r) we obtain 



lim Mir) — lim 



■ u'(r) 



2y/u v 2 



= 0, 



(84) 



(85) 



justifying the assumption M(0) = 0. 

The wave function R diverges logarithmically at the origin. This singularity is weak 
enough that terms proportional to R 2 in the components of the energy-momentum tensor 
(|17fl are integrable. But terms of the form (R') 2 /v diverge as r -4 and are not integrable. 
These nonintegrable terms cancel out of the expression 2T ° — T so the total energy is finite. 

Substituting expressions l|82|) into Eq. (|18|l for T and using R = 87rT gives 



R 



for the dominant behavior of the scalar curvature as r 
diverges at the center of each geon. 



(86) 



0. Thus the scalar curvature 



4 Implications of infinite curvature 

Because the scalar curvature diverges at r = it is not possible to establish a locally flat 
coordinate system there. So the locus r — must be excluded from spacetime, leaving a 
hole in the spacetime manifold. |10| The presence of this spacetime singularity raises the 
following questions: 

1. Do physical parameters (such as mass and charge) diverge? 

2. Do arbitrary boundary conditions arise at the singularity? 

3. Can particles and photons encounter the singularity and, if so, what happens to them? 
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Questions 1 and 2 have already been addressed: the geon mass, charge, angular momen- 
tum and action are all finite, and the wave function is normalizable; and the conditions 
Q 2 (0) = and M(0) = are not arbitrary boundary conditions at the singularity — they 
are consequences of the model [see Eqs. (|51t and IJ85J]. 

The answer to question 3 is not so clear-cut. In a classical (non-quantum) theory one 
considers spacetime pathological if the world line of a freely falling test particle (a point 
particle of negligible mass) does not exist after (or before) a finite interval of proper time. 
Such pathological spacetimes are said to be timelike geodesically incomplete. jlO) In a quan- 
tum theory, however, no real particle can serve as a test particle in the neighborhood of a 
singularity, since such a particle would have to be so small (short wavelength, large energy) 
that it would itself dominate local spacetime structure. Perhaps the concept of geodesic 
completeness can be extended to the quantum realm but at the moment no consensus 
exists on how to do this. 

So we will not be able to answer question 3 definitively. But we will present a classical 
analysis which suggests the singularity is benign. In sections [4.1l and ET2l we calculate timelike 
and null geodesies — the paths of freely falling particles and photons — in the geon spacetime. 
We will find that the geon core is sufficiently repulsive to cloak the singularity and leave the 
spacetime timelike and null geodesically complete. 



4.1 Geodesic equations 

Consider moving along a geodesic with velocity v M . The trajectory obeys 



ds 



2 v 



(87) 



subject to the constraint 



v% = C- 



where £ — 0, 1 for null and timelike geodesies, and s (proper time for a timelike geodesic) 
parameterizes the trajectory. gp„ is independent of t and cj> so equation l|87l) implies v t and 
are constant along the trajectory. Given the spherical symmetry we can (without loss 
of generality) consider trajectories confined to the equatorial plane (6 = ir/2, v# = v e = 0). 
Thus 



( £ .\ 

—vf 



V -i ) 



and 



( £,U \ 

r 



V ey J 



where r = dr/ds, and e and I are real constants. Then equation QSSjl gives 



where 



V 2 = u 



The evolution of <fi is determined by 



V 1 



, e 



<P = d<p/ds = = t/r 



(89) 

(90) 
(91) 

(92) 



and this equation and H9l)[) completely determine the trajectory. Since u, v > equation 
(|90|l can be satisfied only when e 1 > V 2 , and the condition e 2 — V 2 corresponds to a turning 
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point of the radial motion. We will refer to V 2 as the effective potential (for radial motion). 
The effective potentials for timelike and null geodesies in the spacetime surrounding the 
n = geon are plotted in Figs. Hand |H1 



4.2 Geodesies near the singularity 

As r — > the metric approaches 

ds 2 = u dt 2 - v 2 r 2 dr 2 - r 2 dO 2 - r 2 sin 2 6»<# 2 



(93) 



Using the results of section |4~D we will explore the equatorial geodesies of this metric. 
Let ro represent the minimum value of r along a geodesic. By setting e 2 = V 2 we find 



'•() = 



u e 2 



(94) 



The geodesic trajectory is determined by 



dr r I 



(95) 



which can be integrated to give 



±rn In 



— I -1 

r} 



where we have introduced the scaled radial coordinate r* 
of integration so 4> = when r* = 



(96) 



V2 r and chosen the constant 



Geodesies corresponding to Ip^jl are plotted in terms 
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Figure 8: Effective potential for radial motion along null geodesies in the spacetime sur- 
rounding the n = geon. Curves labeled by £ 2 . V 2 — > as r — > oo. 



of the Cartesian coordinates 



x = r cos i 
y = r* sin <; 



(97) 



in Fig. El Test particles are attracted when r* ^> 1 and repelled when r'Cl. In the limit 
of zero angular momentum (£ = Tq = 0) the geodesic is coincident with the positive x axis, 
corresponding to a test particle rebounding directly backwards. 

Now let us calculate the behavior of s as a function of r for particles and photons 
following the geodesies of Fig. [5] From Eq. 190(1 we obtain 



[rfy 



e 2 -V 2 
u v 2 



l-(^ 

r 



(98) 



where 



Integrating l)98|) we obtain 



u v 2 



(99) 




(100) 



where we have introduced the scaled parameter s* — w 2 ks and chosen the constant of 
integration so s* — when r* — 7q. Equation (|ll)0|l is plotted in Fig. ^| for various values 
ofr*. 

Figures and ^3 suggest that the core of the geon behaves like a repulsive potential 
barrier, effectively cloaking the singularity. Freely falling test particles with zero angular 
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Figure 9: Geodesies associated with the metric Ij98|l near the center of a geon. See Eqs. 19711 . 
Curves labeled by . 




Figure 10: Parameter s* versus radius r* for geodesies associated with the metric (|93ll near 
the center of a geon. Curves labeled by 7g. 
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momentum (£ = r$ = 0) do encounter the singularity, but they rebound — their world lines 
do not vanish. The spacetime is timelike and null geodesically complete. 

One may object that any test particle encountering the singularity will be destroyed by 
infinite tidal forces. But keep in mind that no real particle can serve as a test particle at 
the Planck scale, so infinite tidal forces at the singularity do not necessarily pose a physical 
conundrum. 

5 Discussion 

The use of boundary condition to force correspondence between Eqs. (|61() and (|63|l can be 
motivated in the following way. Consider two different localized physical systems A and B 
with zero angular momentum and identical charges and masses. Suppose A is accurately de- 
scribed by the flat spacetime Klein-Gordon equation H59|l while B is described by the curved 
spacetime equation (|22|) . Viewed from a sufficiently large distance each system appears 
as a point particle with a wave function given by (|60|l . Using local measurements of the 
probability density <J"I' an asymptotic observer could determine a and thereby discriminate 
between the particles, despite the fact that both have identical charge, angular momentum 
and mass. It is physically appealing to expect the asymptotic behavior of to be uniquely 
determined by these three parameters. Indeed, this is what we typically mean by the word 
"particle" . By selecting ui to force a to be identical for particles A and B we guarantee that 
the asymptotic behavior of \J , \I r depends only on charge, angular momentum and mass, and 
is independent of internal details. 

At the Planck scale many physicists expect classical notions of spacetime to fail — a 
concept conveyed by the phrase (also coined by Wheeler) "spacetime foam". A successful 
quantum theory of gravity would, it is thought, flesh out the details of spacetime foam and 
erase the singularities associated with point particles. 

The geon model developed here suggests a different perspective. Classical spacetime is 
assumed valid at the Planck scale and point particles are replaced by eignemodes of a quan- 
tum field. Singularities reminiscent of point particles remain, but they do not disturb the 
geodesic completeness of spacetime. Multi-particle systems would, presumably, correspond 
to multiple excitations of the geon field |15l |3| and all particle interactions (including those 
involving wave-packet reduction) would ultimately derive from the action Q). In short, the 
geon perspective replaces the search for a "quantum theory of gravity" with the search for 
a "gravitational theory of quanta" . 

Some questions for future investigation come to mind: 

1. Do charged solutions of Eq. (|69fl exist? 

2. Do solutions with nonzero angular momentum exist? The spherically symmetric trial 
solution [Eqs. (|13|l . (|14|l and (|15|l ] could be replaced by one with axial symmetry. 
However, the trial solution would now involve 9 and <fi as well as r and t, and the 
metric tensor and vector potential would have additional nonzero components, so 
the solution would be much more challenging. Furthermore, it is not obvious what 
functional form should replace metric tensor . 

3. What is the physical significance of the angular frequency oj = Mj\fY r . 

4. What is the physical significance of the stationary value of the action and the action 
per cycle? 
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5. The action is, arguably, the simplest which includes gravitation, electromagnetism 
and quantum mechanics. But there is no compelling reason to regard it as correct. 
p. 144] It would be interesting to see how other terms in the action (such as a 
cosmological constant or conformal coupling P p. 116]) affect geon solutions. 

The geons described here are far too massive to correspond to any known particle. They 
would interact gravitationally with ordinary matter so they appear to be candidates for 
dark matter. A number of workers have considered Planck-mass particles as dark matter 
(see [5] and references therein) but it is not clear whether such models can be successfully 
incorporated into standard cosmology. The density of dark matter within a galactic halo 
is thought to be about 0.3GeVcm -3 pQ so, if all dark matter is composed of Plank-mass 
(1.2 x 10 19 GeV) geons, the local geon number density is about 0.3 x 10 -19 cm -3 . This tiny 
density and weak coupling to ordinary matter would make the detection of such particles 
difficult. 



A Notation and conventions 

Our notation follows Dirac 7 . Greek indices take on the values 0,1,2,3 and repeated 
indices are summed. The spacetime coordinates are x M with x° = t and the metric signature 

is H . 

The curvature tensor is 

R%y/3 = — r Q M „ i( 3 + r a M/ 3,„ — v liv v a a p + v a npT a av , (ioi) 

the Ricci tensor is = R a M „ a and the scalar curvature is R = R M M , where 

is the Christoffel symbol and a comma preceding some lower index [i denotes partial differ- 
entiation with respect to x^. The metric tensor g^ v is symmetric, its contraction is 

g^g MI/ = g/ = 4, (103) 

and, for notational brevity, 

V=^/|det(g MI/ )|. (104) 

A semicolon or colon preceding some lower index /i denotes a covariant derivative with 
respect to x M . A semicolon denotes the covariant derivative of general relativity. For a 
scalar, covariant vector and contravariant vector 

v = v — r Q v 

v% = v^ + r^v (105) 

and the procedure generalizes in the usual way to tensors of any rank and mixture of 
covariant and contravariant indices. A colon denotes the generalized covariant derivative 

U^ = U iM + iQA M U (106) 

for an arbitrary tensor U. 
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B Asymptotic wave function 

In the asymptotic regime [see Eqs. and JSSJ] 

§ = Q 2 r- 1 +0{r- 2 ) 



(107) 



and 



u = 1 - 2Mr~ 1 + Q 2 r~ 2 + 0(r~ 3 ) . 
By expanding u^ 1 in terms of r -1 as r — > oo we obtain 

u = 1 + 2A.fr" 1 + (4Af 2 - Q 2 ) r~ 2 + 0(r~ 3 ) . 



(108) 



(109) 



Consider a trial solution for R of the form 



R = Roo [l + ar- 1 +0{r- 2 )]r 



l-er 




(110) 



where a is a real constant. Substituting the above expressions into Eq. 11' I'D and expanding 
in terms of r _1 as r — > oo gives 



Equation I jllip must be satisfied term-by-term so k = and we obtain Eqs. Ij62|) and (|63|1 . 
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